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Te. INTECLUCTION 

A Wariety of acotstic models exist to predict transnis- 
Sion lcss. Each of these medels contain inherent strengths 
and weaknesses. All have shown poor results in a Sshaliow 
water e€nvircnment due to difficulties at the ocean Ectteon 
where there is an interface ketween media of differert sound 
speeds ard densities. 

Since its introduction (Hardin and Tappert,1973), the 
Farakclic wave esuatico has been a widely accepted uzeans cf 
soluticn fer accustic propagation. The earliest rrecrams 
used 4 split-stef Fecurier transform algorithnz to scive the 
paratclic e€guation (FE). Several other solution tecnonijues 
have reer develored jfrimarily to overcome difficuities that 
cccur sken the Fourier transform encounters an ift¢rtface 
kerween different media (Lee ard Botseas, 1982 and £EcDaniel 
and Lee, 1S€2). 

Arn alt€rnative sclution technigue that uses an igflicit 
fianite-difference (III) alcorithm was developed by lee and 
Fapacakis (1979). {has methed ancorporates aprrogriate 
interface ccnditions and allows solutions in shallow water. 
Starting with the IFT algorithm, Jaeger (1983) developed a 
computer mcdel to fpredict transmission loss and accustic 
Fresstre -tased cn user specified bottom topfography ard a 
Single scund speed frofile. This computer model uses the 
Mathenxatical treatmert of the horizontal and sloping irter- 
faces developed by #&cLbaniel and Lee (1982) and Lee and 
McDaniel (1583). It also utilizes several design treatures 
and mwetkcds inccrporated in an earlier computer [frcgranm 
develczjed Ly Lee and Eoteas (1682), and a PE computer model 
geweleped Ey Breck {1578). 
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Tbe JED program freserves continuity of pressure and 
contixruity cf the ncrmal ccmgcnent of particle velccity at 
an interface between media having different scund speeds and 
densities. This feature makes the prcgram uaconditicnally 
stable and better atle to handle tne bottom Ecundary 
conditicre. 

Since its deévelctment, the IFD program has net ceen 
Kigorcusly tested. This thesis analyzes the jfrograg's 
Ferfcrmarce in an idé€alized shallow water environment. The 
envircnuert includes a simple sioping sand fottom Feneath ar 
isospeed water field. The analysis pdegins by comparing the 
ITPps ELI¢Qleticns ae predictions from Jenrscn and 
Kuperman's (1980) FE model and Coppens, Humphries and 
Sander's (1684) imace model. The analysis also includes a 
comperiscn of the acdel's estimated transmissicn icss 
contcurs with expectations based On sitple Fhysaical 
reascninsg. Finaliy, tke thesis describes an attempt tc mcdel 
a shallcw sloping bottom in an experimentai tank. Mhe tank 
contains a sand bottcm sculptured with a ten degree silo;ye. 
laboratcry weasurements ci tne pressure field are taken at a 
frequency cf 100 kHz for ccmparison to the predicted ;res- 
sure field cenerated by the IFD. These compariscns cf the 
IfD predictions with cther acdel eStimates, theory, and 
laboratory measuremerts jgiye an indication of the iFI's 


pericrmance in a Shallcw water environment. 
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Aw. EACKCECUND 


Ac iggydicit finite-difference solution technigue te the 
Fararclic eéyuation has been studied and refined ty many 
autnoors. The nistory of this development is explained in 
detaii by Jaeger (1982), but merits review in order te gain 
a perspective on an aralysis of the IFD computer model. ine 
Farakclic equation is an approximation to the ellirgtical 
wave e€guaticn. The first means of solving the PE used a 
split-step fast Fouriet transform method as develcped by 
Taeeert and Hardin (1973). This method reguires feéeriodsic 
rFoundary conditicns it depth recause of the finite Fourier 
Meeamertorm and handles this constraint cy introducain¢ ar 
artificial horizontal fressure release bottom belcw tae 
actual physical bottcn. This method of implementing an 
artificial bottom was inccrferated into the earliest PE 
models developed by Jénson and Krol (1975) and Brock (1578). 

ELrcrs in this split-step Fcurier transferm methcd were 
found tc be progporticnal tc the horizontal range stefy and 
the second derivative of the index of refraction. The seccnd 
derivative of the irdex of refraction tends to te laige 
across tke ocean bottcm interface. Another problem with the 
split-step method is that it does not consider density 
Gifferences between two difierent media at an interface, 
which ixfiluences the reflection coerficient. For these 
reascns the split-ster fourier transform methcd preved to be 
Fooriy suited for a stallow water environment. 

The IFC soluticr method was introduced by Lee and 
EFapadakis (1979) as ar alternative to the split-sterp methcd. 


The IFD method emplcys a seccnd order central difference 


/ 


formula to solve the FE in the £crm of a@ tridiagonal fatrix 
Althcrgh the first version of the ifD did handie discontiru- 
ities in the scund speed profile, at didynot consider 
density discontinuities. In 19&2 HcDaniel and Lee intrcecuced 
a meétikcd tc harndie a horizcntai interface of ditferent 
densities. in’ 9S 3; they extended tneir treatzent to 
include a sloping interface. It is this version of the IFD 


that is used in Jaeger's computer proyran. 


E. TEE CCHEUTER MODEL 


The JFL computer program consistS of 42 Main progra ee 
twenty sukrcutines. tThe precgar utilizes a modular ccnstruc- 
tion sc that each of the various subroutines are called fron 
the main frogram to ccmpiete a Specific calculation cr furc- 
tion when trequired. The TFL iS run interactively fren 
user generated PD put ea ie thace Genta. values for 
freguency, one souraq speed profile, a bottd)&@ pro iagier 
sourceyreceiver depths, attenuation coefficients fer rctn 
the water and the bottcem, and several other input parameters 
that tell tre the precram where to obtain a solution within 
the field. The procram initiates the calculations assuming 
an initial Gaussian fressure field and an artificial [;res- 
sure release surface at a user specified deptn. 

Atteruation in tkoth the water and tne sedinrert is 
handled as complex indices cf refraction. An artiiiciai 
attentation layer is estaklished beneath ture sedinert to 
introcuce attenuatior abcve the artificial pressure release 
surface. The actual fagnitude ox this enhanced attenuaticn 
is calculated using ar eguation derived by Brock (1978) for 
use in his FE computer nodel. 

lke IFD pregram steps along the vspeetilege oot vou, per 
and gakes calculaticns dcwn tarouyh the water/sediment 


column at €ach user sfzecified horizonal range. The rpicgran 


/ 


requires that tne Ecttcm intersect exactly at a vertical 
Simi pcirt. Po dee Coupe tO Los lO PMN hone Ol Lhe pretax 
autcmatically calculates the range step to fulfill this 
reguiremzent. This ccmputer generated range step is always 
idess than cr equal tc the user provided range step. As the 
Slope ci the bottom increases, the range step must décrease¢ 
and gore calculations and computer time are reguirec to 
solve tke entire pressure field. In some Situaticns an 
actual mcdification cf the user inputted bottom is reguired. 
Thais cccurs with a very gently siopiny bcttom when the 
reguired rarge step exceeds the user specified rancse step. 
Here, the program automatically modeis the bcttcm as a 
series cf level and sloping sections in order to ensure the 
user cené€xrated range step is net exceeded. This mcdification 
cf the Ecttcm is always less than or egual to one-half the 
vertical grid spacing and the model issues a warnirg tc the 
user cf the mnodificaticn. 

Ecth printed and graphical output are provided fry the 
IFD. The fErinted cutput prevides transmisSicn loss and the 
real and iméginary ccufonents cf the pressure field at eack 
depth fcr a specified hcrizcrtal range. The graphical cutjut 
is a glict cf transgission loss versus range at the user 


specified receiver defth. 


C. #ECDEL EROBLEMS/MCLIFICATICNS 


As this study of the IFD fregram progressed, it became 
necessary tc modify certain aspects of the model. Mcst of 
these mucdifications sere necessary to alter the fprcgiran 
cutput irtc a mcre desirable fcrno, put a tew were impie- 
mented tc correct picgramming deficiencies. Although this 
secticn cf the thesis discusses the Earliest model rurs, 
these results are rect presented in detail, ‘Gig  feedkyy 


discussed in general terms recause they were obtained Lefcre 
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the ccuputer model sas fully mccified. All modificaticnrs 
here gadcde cnly aftezr careful analysis of multiple model 
runs. It is impertant to realize that all results generated 
for ccmparison to otter models and laboratory measurements 
came frcnr a fuliy mocified version of the IFD. 

Since the ultigate okjective of this study was to 
compare 1FD model results with laboratory measurements the 
Frogram@ was Eixrst CUI with irput parameters waren wexvae ws 
modeled ccnditicrs in tne tank. The experimental set uf, 
explained later in great detail, consisted of a tank tkat is 
approximately twe meters in lergth, one meter in depth, with 
aten cecree slcping sand Ecttcm and a maximum water derth 
cf 35 centinzeters(cm). Based on a test case run ky Jaecer 
(1983) cf a simple sloping bctton, and by results shcun by 
Jenscr ard Kupergan (1980) for propagatich in a wedge-shajyed 
ccean, it was expected that there wouid be certain recccniz- 
able patterns in the predicted propagation patterns. 
Specificaliy, since tke speed cf sound in the bottom exceeds 
tnat in the water {fast bottcn), tke simple slofing [Ectton 
supperts trapped acrmal =acde propagation (Cofrens and 
Sanders, 1981). As the acoustic energy traveis upsicpe 
toward tke apex, successively lower modes are cut cfi anda 
the e€nercy contained in these acdes is transmitted intc the 
bottcno. The range from the apex at which energy cf the 
lowest mcde is transritted inte the bottom is referred to as 
the dump distance anc is a function of waveiength, wedge 
angie, and the ratic cf sound speed in the water tec the 
sound speed in the sediment. An empirical equaticn tkat 
defines this dump distance was derived by Copfrens, Sanders, 
Toanncu, ard Kawamura (1978); and was used to identify the 
expected ranges of these dump distances for the civen 
scenario déscritked aicve. 

Ite initial unmodified IFD run used the parameters taken 


from the tank and showed fo recognizable patterns ir the 


eeguetic f1¢1d. There was no ckservakle decrease in trauns- 
fissicn lces at the various dump distances as expected. 
Father, results indicated widely fluctuating patterns in the 
acoustic tield that appeared inconsistent with both rrevicus 
studies and sizple physical reasoning. Upon closer analysis, 
it was disccvered that althcugh the prograz was designed to 
re incependent of scale there are several logic statezents 
that are net itplemerted if the user provided range ste;z is 
jess thar cne meter. Because the loyic statements aren't 
Satisfied, the NEWSEG and NEWMAT subroutines (Jaeger, 1963} 
ace nct called ccrrecitly. The NEWSEG subroutine initializes 
the Ecttcm siope anc the NEWMAT subroutine computes matrix 
Glomernts for the precram. Obviously, errors in these two 
Frograg functions seriacusly distort results. Because of this 
systematic error in tke program it became necessary tc scale 
up ail tank parameters. Ali distances were scaled uf ty a 
factcr of 1000, and frequency was scaled down by a factcr of 
1000. Careful analysis reveals that all input farameters are 
afurncticn cz either distance cr freguency, so this scaling 
Froduces results that model these expected in the tank. 

The second modification cf the IFD was reguire¢c to 
provide a three dimersional graphics display. As discussed 
€arlier, the IFD prevides a transmission loss flot versus 
Tange ata single cepth. Hcwever, to study the saodei 
predicticns in greater detail it was felt that a twe dimen- 
Sional aralysis cf tke model estimates would be mcre tean- 
ingfuvl. As a result a transmission loss contouring frcgran 
was developed. The program {Appendix B) displays transmis- 
Sion lcss contours for range versus depth. Use cf the 
Gentcur ;lct requires that transmission loss vaiues gener- 
ated ky tke IFD be sent to a data disk used by the ccrtcur 
Foumiure. Ic facilitate this transfer a dummy variakle (TifD) 
was estaltlished in the PRINI2 subroutine to store the trans- 
Dissicn loss valves and then these values are written tc the 
Gata disk at the end cf the main frogran. 
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Anctker modification of the IFD output was reguired to 
change tke reai and imaginary components of the [fressure to 
a Sincle pressure auylitude tagnitude. It was felt that 
Gealincg with the pressure magnitude was easier ard mcre 
gGeanincful than witk the real and imaginary comecnents of 
the ,~ressure field. This ccnversion was done in the PRINT2 
Ssubroutire and established a new variable (PRMAG) tc repre- 
sent tke pressure macritude. 

A final modification in tke computer prcegram was mace 
due tc a Suspected errcr in the computation of the attenua- 
tion in the artificial layer. FPhysicai reasoning dictates 
that proper impliemertation of the artificial attentuation 
would result in a steady drcp off in acoustic pressure with 
deptk thrcoushout the artificial layer, With pressure drop— 
Fing tc zero at the pressure release surface. IPL Kedcet 
results cn the other kand shcwed wide fluctuations in fres- 
sure sith depth in tke layer and then only a minizal fall 
cff at tke fressure release surface. The egquaticr in the 
IFD tbat actually computes the magnitude of the attenuation 
in tke artificial layer was taken directly from Exrcck's 
{197&) FE mcdel (Jaecer, 1583). However, Brock's equation 
WwaS derived with feet as the unit of measurement r¥hile 
Jaeger's acdel is derived witb meters as the unit cf treas- 
urement. whith this 1n mind, Jaeyer's eguation skEcula Ee 
apprcximately a factcr of three larger than Brock*s equatacn: 
to ccrrect for the difference in units. To correct fcr this 
error tke ecuaticn tc calculate attenuation (ATT(I)) in the 
NEWMAT sukrcutine was increased by a factcr of three. Wen 
this ccrrection was implemented, the large fluctuaticrs in 
Fresstre with derth were eliminated. The expected dror cff 
in ~ressure with deptb and the fall off of pressure tc zero 
at the pressure release surface were noted. 

A listirg of the revised IFD computer program witk all 
Godificaticrs can be seen in Appendix A. 
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Te. MCCEL VERIFICATICS 


lowmrig¢e: Of Whe modifications to the LFD cerputer 
Frogragz just discussed, it was necessary to ensure that the 
the chances thenaselves did not introduce errors irtc tke 
model. ‘So as a first step the modified IFD frogram was run 
for cre of the test cases used by Jaeger in his criginai 
work. Ibis case analyzes propagation in an environment that 
moves £xrcm deep to skallcw water. This particuiar envirca- 
ment is depicted in Figure Z.1 and was chosen because it was 
very Sinilar to the simple sloping bottom in the tank exjper- 
iment. A scluticn is cktained for a bottom with an ufsicpe 
cf 8.5 cegrees. An isospeed water field is used wati scund 
speed set at 1500 oys. Scurce freguency is 25 Hz. ihe 
scenatic has a maximunz depth of 350 meters anda rancge cf 4 
kilcneters. Both the source and the receiver are set at 25 
Beters. 

The results using the rodified IFDP program were the 
same as Jaecer's for this test case above the artificial 
attertuaticn ilayer. As explained earlier, the prccrai was 
modified tc reflect higher attenuation in the artificial 
layer, <s0 as to preyerly reflect the effect of attenuation 
in this regicn. The modified program results did ret show 
the larse fluctuaticrs in pressure with depth in the artiti- 
Cial layer but ratker the gradual decline in [fresstre 
towards a value of zero at the pressure release surface. Eut 
above tke artificial layer the modified program results were 
identical with those achieved by Jaeger with the crisginal 
Frogran. Apparently, the miner caanges in the fregras 
Gesigned tc improve cn the form of the program outfut dces 
bot binder the mcdel's ability to achieve a soluticr in the 


upper sé€dizent and water Column. 


1€ 


FREQUENCY —- 2.0 HZ 
SOURCE DEPTH —- 25 M. 


SEDIMENT 


ARTIFICIAL LAYER 


RANGE(km) 





Figure 2.1 Jaecer's Deerp-to-Shallow Water Case. 


uperman Model Run 





The seccnod attempt at verifing the IFL- macdel 
lnvolved a compariscn of mcdel] results with those achieved 
with a FE mcdel designed by Jensen and Kuperman (1980). The 
Jensecr and Kuperman Mcdel (JKM) uses a split-ster solution 
techrigue. Compariscr with this particular model was chosen 
kecause it is one of the few that obtains a soluticr ir two 
dimensicns. Most acorstic mcdels obtain a solution at only a 
Single depth. In addition, Jensen and Kuperman made their 
model rurs in a simple sloping ccéan bottom envircrment very 
similar to the scenario of interest modeied in the tark. 
This €rvircnment is depicted in Figure 2.2 and features a 


gently slogping bcttcn of 2.2 degrees. The water column has a 


1S 


Uitemera Speed=ot 1500 gmys. The source is placed just kelow 
the midpcint in the channel at 112 meters and is driven ata 
frequency cf 25 Hz. The maximun depth in this scerario is 
z00 meters and has a maximunx range of 12.5 kilometers. 


The Jensen ard Kupérman results were expressed as 


EREQUENC Yee 2.5 HZ. 


SOURCE DEPTH -112 M. 
0 


WATER 
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Figure 2.2 Jensen and Kuperman Sloping Bottom Case. 


transnission loss ccrtours fcr range versus depth. their 
study ccncentrated cn transtission loss patterrs in the 
sedinzent, Eut results were cbdtained both in the sedizert and 
in the water. These results are compared to the IFD results 
meee OP. Spo. 0, 1-5, 10-0, and 12-5 Kijiometers. iInese 
ranges were chosen for analysis because the JKM results 


showed the greatest variaticn in transmission loss with 


ZN 


depth and thus, makes for a more meaningful compariscn with 
the JIC results. 

The IFD and JKM results can be seen in Figure 2.3 
through Figure 2.7. In all the figures, the IFD €stizates 
are Sicwrk asa solid curve while the Jensen and kKuperiman 
results are depicted as circular points. MThe first analysis 
is at a rance of 2.5 kilometers. From Figure 2.2, it can be 
seen that tke depth at this range is 200 meters and ig ira 
flat kottom regicn. From the results shown in Figure 2.3, 
it is ckvicus that [Ecth mcdels obtained almost identicai 
results frem the sturface dcwn to a depth of abcut 300 
Meters. There are differences between the two sets of 
Fredicticns from the ccean Eottom to a depth of 100 meters 
kelow this foint. Eelow 300 reters the Jensen and Kufeérman 
results shca a very slight increase in transmissicn icss 
(IL) watk depth. Tke IFD alse Shows an overall increase in 
transgissicr loss with depta rut with several fluctyaticns 
in txrarnsenission loss and a marked peak at about 300 metetLs. 
So in general, the results frcm the two models have the Sate 
generai terdencies altkough the IrD appears to shcw greater 
detail in results near the water/sediment boundary. 

Figure 2.4 stews results at 5.0 km in range. Here, 
the water depth is still 300 meters and marks the very 
feginrirg cf the sloping bottcm section. Fer this rance the 
JKM predictions are criy available to a depth of 300 meters. 
Tae results are nearly identical with those cbtained Fy the 
Iie Ecth models show a relative minimum in transmission 
loss at a depth of 1€0 meters and then a gentle increase in 
TL with depth. 

Tne model results at a range of 7.5 km are seen in 
Figure 2.5. At this range the bottom depth 1s akcutwisadg 


meters anc the bottcm is slcping. Ine medels show the 
greatest difference at this range. Both models fpErcduce 
nearly identical results in the water column, but beneath 
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Pigure 2.4 IFD and JKM Conparison at a Range of 5.0 Kn. 
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Figure z.5 IFC and JKM Ccmparison at a Range of 7.5 Ka. 
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150 meters in depth the mcdel estimates ceyin tc slew 
éifferences. The Jensen and Kupercman results Show a4 gradual 
increase in TL with depth below the ocean Lottoa. ide wep 
results cn the other tand, shcw a gradual decrease in tLrans- 
missicn lcses frem 1£&C meters to 300 meters, with a relative 
minigun at 300 meters and then an iacrease in M@L frereath 
this de;th. Althouck there are some differences in srodel 
estinates rear the [cundary in the flat pottom secticns 
previcusly discussed, the differences appear to ke greater 
at the kcurdary in this slcfing bottom case. Ir thcse 
regicns away frog the Ecundary in either the water cclutmn or 
in tke se€dineat, hcweveLr, Ecth models produce sigilar 
results. The JkKM's difficulty in obtaining an accurate sclu- 
tion at tke sedigent [cundary is not totaily unexpected. ihe 
model tsé€s a variaticn of the straight split-step sclution 
technigué Simiiar tc Brock's computer model (Jenser and 
Kupergan, 1980), trat has characteristically been unsuc- 
cessful in cbtaining @ reliable sclution aeéar a boundary. 
The resuits at a range of 10.0 km are seen in Figtre 
Or At this range the water depth 1S apfroxitately 100 
geters and the bottom is again 1n a Sloping region. At this 
Tange Jensen and Ktperman results are only available to 
about 15C meters in depth. For the data availatle, the 
models predtce nearly identical results. Beth models fredict 
a minintug at about £C meters in depth and then an almcst 
linear increase in IL with depth. The IFD results also 
reflect @¢ transmissicr 1lOSS faximum at about 200 meters, and 
then a slight decrease of IL beneath this depth. Pa ta eat 
these derths are not availatle from the Jensen and Kuferman 
Eun In dight of tke results cbtained at 7.5 kt £Or ea 
Sslopirg rettom case cne might expect that the results fcr 
the two acdeis wculd show differences near the ocean Lottcn. 
However, Since there is cnly one Jensen and Kupferman 
Fredicticn available rear the rtottom for this range, it is 


difficult tc nake any definite conclusicns regarding differ- 
ences it gcdel perfcrmance at the boundary. 

Fesults at tke apex (range equal to 12.5 kz) can Le 
seen ip Ficure 2.7. Since all estimates are nade in the 
Sediment at tnis rarge, Prmelo 1S00 Doundary tc GcLtend 
with. At this range the twce models snow the test acsreezert. 
Eoth models show a gradual decrease in transmission licss to 
aderfth cf about 300 meters and then a gradual increase of 
pieeetl derth. 

In general, there appears to te good acreenent 
Fetweex tke IFD and the JKM godel results in regicrs act 
iantluenced Fy a bouncary. in koth the water and the deep 
sediment the twe models preduce simiiar results. Aitkougn 
this is ret conclusive evidence, these sinilarities sucgest 
that the ITF2 can successfully model accustic propagaticn ir 
these reciccs. Near tke waterysediment boundary however, the 
Jensen and Kuperman ard IFD predictions show marked differ- 
ences. These differences afpear to intensify as the Ectton 
reccnes mcre sophisticated. in general, the Jenser and 
Kupergzan results dcret sé€enm to snow the detail the IFD 
results do. Considering the different solution technigves 
emplcyed Ey the two acdeis, tnese differences in resuits are 
expected. 
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The third attempt at verifying the IFD Lperfcriance 
was dcre Ey comparing results with an image tueory rocei 
derived ty Coppens, Humphries, and Sanders (1984). This 
Coppens, Humphries, and Sanders Model (CHSM) uses a saddie 
Foint appreximation to an image model tc soive fcr the 
deguetic tieid. Beth prograzs were fun for the scaled 
ecenaric mcceled in the tank €xperiment. This scenario is 


depicted in Figure z.8 and features a ten degree slcping 


Ze 


kottcz. The maximum water déefth for the run is 350 meters 
and the raximum range is twe kilometers. The source gener- 


ates a 100 EZ signal and 1S Set at a depth vor sence cc- 


FREQUENCY -100 HZ. 
SOURCE DEPTH -175 M. 


WATER 


SEDIMENT 


ARTIFICIAL LAYER 


he 
RANG=(km) 





Figure 2.8 10 LTegree Sirple Sloping Bottom Case. 


The image model fFrovides soluticns only at the apex. 
The ccnrpariscn of fredicticns at the apex(range egual tc 2 
kilometers), can be seen in Figure 2.9. In the Figure, the 
IFD values are flotted as the solid curve while the image 
model values are shown aS Circular points. The predcicticns 
are Sskcwr as values cf normalized pressure amplitude for a 
Given depth. The different models produce pressure anplitude 
values in different writs, and thus, had to be normalized to 


wake a CcC@fearison of values possible. This normalization was 


cone fcr each model ty dividiny the presssure anplittée by 
the maxigua ampiitude predicted by the modei. 

ihe results suggest that even though the tcdeis 
Fredict sinilar larce scale trends in pressure aunplitude 
with depth, there are several differences in détail. The 
Flot shcws that both the IFE and the igzage medGel predict az 
almost linéar increase in amplitude with depth Corn tea 
specific maximum and then a Slow decrease in aupzlittde 
Feneatrk this maximum. The IFD however, shows tae faximum at 
about 15.5 meters in depth while the imaye model maxinun is 
deeper at 1 meters. Beneath this maximum the IFD shecwe a 
fucn sStarper decreasé in amplitude than the image rcdel. the 
trend in the izage mcéel data appears smooth, wnile the IFD 
curve reflects several small scale fluctuations. 

It appears ttkat both models predict sSinilar large 
scale trends in [resstre amplitude for this scenario. It is 
difficult, if not imgessible, to account for the differences 
ain détail. AS a Minigun, hcwever, this comparison indicates 
that the two sets cf predicticns are consistent with cne 
another and can be considered reasonable in this sShailow 
water envircnment. 


4. Cctpariscn with Physical keasoning 


Frex modél cezparisons it appears that the ifr at 
least makes a reasonatle estimation of acoustic fields ira 
Shallcw water enviroruyent. IFD results are also analyzed in 
CCMyparisccen with basic physical reasoning and theory és a 
fourth atteapt at model verification. This analysis ¢xarines 
IFD transmission loss centcurs for the simple ten degree 
sloping cceéan scenaric seen in Figure 2.8. The scenario is 
two kilcuzeters in rance with &@ maximum depth of 350 meters. 
The tke Source generates a Signai at 100 Hz and is placed at 
175 meters in derth. 
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Figure 2.9 IFC and CESM Comparison at the Afex. 


Itarsmissicn loess ccntour plots are disrplayed in 
fguce 2.10 througk Figure 2.13. Vv ioe CONCOURS are 
Peepeecsed ai decibels(db) and are shewn as a functicn of 
depth versus range. In ali figures, the ten degree slczing 
rottcna is shown aS a solid unlabeled line. The contcur 
Fiots are cisplayed over fcur range subsections due to 
computer <sraphics lanitaticns and to emphasize different 
Significant features ir the field. 

jhe first ficure shows transmission ioss ccntowurs 
Irom the sctrce to a range of 600 meters. MACBEECHtCurLrS we 
the first 250 meters appear very symmetric, increasins 
cutward fron the sotrce in a fattern that resembles sfher- 
ical spreading. Since the IF2 program assumes a Gaussian 
starting field this c«érly pattern is expected. 

Frce@g@ a range of abcut 300 meters to 60) meters the 
field in the water appears tc be dominated by a surface 
reflecticn jfattern. From the spotty appearance cf the TL 
ccntcurs in this regicn there 1S an indicaticn of ar irter- 
acticn cf surface reflecticn and bottom reflecticn cn the 
contcurs. It is fossitle to ccmpare transmission lcss 1axina 
with rcedes in the surface interference pattern. Eased on 
surface interference theory these nodes should cccur where 


(Kinsler, Frey, Coppers, and Sanders, 1982): 
SIN(khdyr) = 0 
cr in cthexr words; 
hd we = yay ti 


where; 
=a0, 172. see 
= Range 
= Wavenumber(2'y,) ) 
Scurce Depth 
= Depth cf Node 
= Wavelength. 
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Figure 2.10 
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IFD TL Contours(db) frogs 1350 to 2100 Meters. 
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Ey maripulating the akove €iudation 2t is weec cr eet ce ce 
for tre depth at which the necdes should be observed. for 
this particular scenario at a range of 604 meters, rodes 
should cccur at integer multiples of 25 meters in Cefth 
(25n). Figure 2.1C shows that these transmissicn I1css 
maxima dc exist as expected, every 25 meters in depth at the 
statec range. 

Eased on simple physical reasoning one weuld excrect 
refractive tending aleng the ocean bottom due to density 
differences between tke water and the sediment. The ccntcur 
flot reflectS acharcgée in FEasic pattern at the interface. 
There is a kending of the TL ccntours that suggest a refrac- 
tive infiluerce. 

Figure 2.11 shows transmission loss contours from a 
range of 60€ meters tc 1350 meters. Again in this recicr the 
water aprears to be dcminated Ey surface reflection. Sclving 
for the ncedes in this surface interference pattern at a 
range of 1350 meters, it is found that taese nodes stculd 
appear every 55 meters in depth. From the figure it iS again 
seen that Tl maximums do occur every 55 meters in depth as 
anticirated. As in the first figure there 1s a change in 
the kasic TL pattern at the ocean tottom. The tLerding 
appears fore accentuated thar in the previous figure, fut 
still suggests the influence of refraction at the cceéan 
Forreenu. 

Figure z.12 displays transmission loss frem 1250 
meters tc fast the apex at a range of 2100 meters. In this 
figure, tke dominance cf the surface reflection mechanism is 
less ckvious and the TL patterns become more complicated. It 
is it this region that the influence of trapfed nceraal acde 
Fropacaticn can te seen. As discussed earlier, as accustic 
energy travels up tte slope tceward the apex, noruzal modes 
are cut off and energy is transmitted into the Ecttcz. 
According te adiabatic normal mode theory, mcdal separation 


iS "farce deyendéent (Graves, Nacel, Uberall, and Zaur, 175 
ememece ;chb= and Sancers, 1990). The gcanye from the apex at 
which tte lcwest mcde is transmitted into the bottom, can be 
calculated using the fcllowing equation (Coprens, Sanders, 
Icanncu, and Kawamura, 1978): 


xX = A/4 SING TANG 


wtere; 
% = DumE Distance Of The Lowest Mcde 
)) = Navelergth 
&, 
G 


Craseical Angie 


Wedge Angle. 


Accoréing tc adiabatic norzal mode theory (Kinsler, frey, 
Coppers, and Sanders, 1982), in deep water (near the scurce) 
the lcwer normal modes are far above cutoff and the adia- 
Fatic ¢€icenfunctions consist ci an integer number cf half 
Sine tnaves with zerc fressure at both the top and rctton 
surfaces. At the cutcfi£ of each mode, the pressure at the 
kottcna mest be maximized, resulting in an adiabatic eigen- 
functicn that contains 1/4, 3V4, and 5/4 wavelengths at the 
respective cutotf distances cf 1X, 3X, and 5X for the three 
lowest ucdes. From Figure 2.14 it can be seen tnat as the 
ncecmal mcdes travel up the sicping bottom, succesive rodes 
ace ferced into the Ecttom at distances where a farticular 
mode reaches a deptk at which it can not longer frecfacate. 
Also frcem the figure it is cbvious that a source set at 
Baid-derth can nct eéxcite the second node. Since thkis 
Farticular ceometry is present in the tank scenaric it is 
expected that the energy asscciated with this seccne mcde 
snould ret be seer. Based on this line of fhysicai 
Teascrinoc, modes shculd be dumped into tne sedirent at the 
first dump distance, TPreoeolipmeear Stanec, Menta od ume 


distarce and so cn. Jf these modes are dumped as descrited 
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Figure z.14 Normal Bode Profogation in a Wedge Shared Ocean. 


then they should be seen in the contour flot as reducticns 
cf trarsgission loss in the bottom at ranges cf 1952 zetels, 
1760 geters, 156€ meters, ard so on. 

Frenr Figure 2z.12 it can be seen that there is an 
cbhvicus decrease in tIansmissicn loss at approximately 1650 
meters (dump distance #1) and 1760 meters (dump distarce 
#5). Ibkere is also a less clearly defined reducticn in TL 
along tke Eottom at 1565 meters (dump distance #9). 

ihe final cortcur plot (Figure 2 oie cnovws trannies 
Sioa lcss contours from 1550 meters to 2300 meters. This 
figure is ar extensicr of Figure 2.12, intended to emphasize 
how clearly the beam at the first dump distance is defined 
in; thes lcot=-2 Ine kéam at the rifth dump distance is also 


visifle rut is not nearly as well defined. Frem this figure, 
eucmecen dlcc See an indication of “a very narrow Lear ir tke 
sedigent at about 18€C€ meters. This distance corresrerds to 
the tkird dump distance (second normal mode). Based cr adia- 
Fatic mwcde theory tke seccrd node is not expected te te 
excited. Hcewever, adiabatic medée theory is only an aprroxi- 
maticr cf rermal mode behavior. This approximation cf nersal 
mode kEekavior keccaes less exact as the bottcenm sicre 
increases and the clcser the source is to the apex. Ine 
appearance cf a narrcw beam at the third dump distance indi- 
cates that there iS a strceng fossibility that the seccnd 
cde is present and tkat the adiabatic approxinaticr is ret 
exact with a ten degree bottcn slope. 

The basic features of the IFD contour ficts are 
consistent with botk physical reasoning and theory. Easic 
surface reflection and kEottom refraction occur where 
expected ard behave as anticifated. Pn thes (Lag, itaeld, 
trapped norgal mcde fropagaticn is observed and can te veri- 
fied with simple dumy distance calculations. The locaticn of 
keams dumped into tke bottcem appear consistent with Lasic 
mode thecry. In short, the transmission loss contours indi- 
cate that the IFD 1s making reasonable predictions cf the 


acoustic field ina stallow water environment. 


See Verification Summary 


RGeleraittacult to say how exact the IFD precicticns 
are fcr a shallcw water envirenment based on these single 
verification technigues. As a minimum it can be said tzat 
the ucdel results are at least consistent with other srodel 
pFredicticnse and expectations tased on Simple fhysical 
Teascnhincg. Model estimates are virtually the same as the 
Jenser and Kupergwan FE model in regions not influenced cy a 
wateryseGizent touncary. Clcese to the boundary the iFD 


results appear to skcw greater detail and variaticn ttan 


4C 


this EE medel. The IFD results also appear consistent with 
the general trends in j;ressure amplitude predictecd Ey the 
Coppéens, Humphries, and Sanders image model. Again dizrer- 
ences were roted in the small scale structure. Finally, the 
IFD Tl cecntcurs verify well with rtasic expectations cased on 
Fhysical reasoning and theory. Surface reflecticn and 
krottcnr refraction patterns are observed as anticipated. far 
field neraal mode frepagaticn can te verified in the ficts 
uSing simuple dump distance calculations. In short, ail veri- 
ficaticr methods attempted, faii to uncover any incorsis- 


tency in IFI performance in a Shallow water environment. 
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A. EACKCECOND 


The waltcr attempt at appraising IFD performance inivcived 
comparirg meodel results with lakoratory measurements. MThe 
shaiicw water envirorment mcdeled in tae tank is very ideal- 
gazed; aten degree sloping sand bottom with an iscvelccity 
watec ccluaon. Althcugh this Scenario appears exterenly 
Simplistic, it is ore that carn ke reasonably mocelec ar the 
mie@ratcry and still approximate conditiors in an actual 
shallcw pater oceah ervironient. The sretnods used tc nodei 
and seasure the accustic field are reiatively untested. 
Indeec, this attempt at Llakoratory modeiing was yperfcrzread 
Got cniy tc verify Ifl predictions, but also to see if the 


envirerment could be successfully modeled in the latkcratcry. 


EF. EXEEFINENTAL DESICN 


lie The Tank 


A fiberglass coated wocden tank was used. The tank 
is 304 centimeters ir iength, 117 centimeters wide ard 95 
centineters deep. Sand filling the bottom of the tank was 
shaped tc foro tke ten degree slcpiag bottog, and 
measurements were taken over a range of two meters. Maxizun 
weer derth in the tank was 35 centimeters, A 100 kHz 
source was placed at mid-channel depth (17.5 centimeters). 


The laycut cf the tark is depicted in Figure 3.1. 


"da 39S yuer Tejyuemtaisdxqg 4r-¢ arnbry 


[nL 0) 


—— — O20 0CG -— ae a ae aaa ta uo cS 4 





(NOLLOE LSVAJGNVS 





oe == = —= —_— = — —_—_— =e = = oe 


HALWAA 3JOuNOS 


' 
J 
] 
) 
‘ 
‘ 
’ 
! 
! 
! 
4 
§ 
) 
5 
, 
t 
' 
j 
( 
t 
t 


IV 





| pe me ee em eH Hee 


. 
OV 
O 
a 


- =. 


43 


A slope of ten degrees was selected fer several 
reascns. ic begin with, even though a ten degree slcfe is 
greater than most ocean bottcn slopes, it is still enall 
enough tc ke considered realistic. Pernaps nost Deon taut, 
the ten degree slcfe was selected because given the 
frequency limitations (discussed later in this chapfter) and 
range jdinitations, this wedge anyle allowed the source tec be 
Flaced many (471.6) dump distances fror the apex. A large 
Dumber cit dump distarces was necessary to simulate a distant 
source. 

The bottcm material used in the experiment was #30 
fine crade sand. The grain size ranged froa 0.15 millimeters 
to 0.70 nilimeters. The sand was treated wita a technivue 
used ty Faexk (1984) tc remove air from tne sediment. This 
technigue used a high speed jet to agitate the sandyxater 
Mixture tc remove tke bubbles and then allowed the sand to 
settle fcr several days before the experiment was initiated. 

Fresh(taf) water served as the medium in the tack. 
To remove air burbles, the water was allowed to settle ina 
settling tank before being transferred to the experimentai 
tank fer use. The water in the tank was periodically treated 
miter Chicripe bleach to prevent the growth o£ bicicgical 
Baterial. 


oe 


ze Signal Generating/Receiving Equipment 


The acoustic signal used for the measurements was 
produced Ey a functicn generatcer, sent through an amplifier 
and then transmitted into the water by a directional trars- 
ducer resonéent at 100 kHz. The dimensions of the active face 
cf£ the transducer were 7.0 cm in width and 2.0 cm in height. 
These dirensions resuited in an approximate beamwidth (ancie 
from the acoustic axis to the first theoretical null) of 
fgeeemeecgecrees in the horizontal and 46.3 degrees ir the 


vertical. The narrcw horizcnotal beam minimized reflecticns 
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from tke sidewalls ci the tawk, while theme verermea! seat 
allowed ccuplete ersonifticaticn cf Jrhnewechanre! pred 
vertical dizensicn. 

Tke signal was received Ly an LC-5 omni-directioral 
hydrcpheré, sent thicughn anh aipitricre apa eitves andere] 
displayed cn an oscillosccpe ang voitmeter. A schematic 
chowing the electronic setup is shown in Figure 3.2. 

A freguency cf 100 kH2 was selected for thc rrac- 
tical reascns. First, to avoid  particie Scaeseering Beate 
sediment, the acoustic wavelength wust be at ieast three 
times largexr than tte grain size of the sedinent (Andeérsor 
and Liektérmann, 196€&). The itargest grain size in the Sand 
was C.C7? cn, so that a wavesentn of 1.45 cm(100 kEz) was 
sufficiertliy large erough tc be immune to this eéffect. 
Second, the 100 kHz freguency and tae properties or the sana 
Frovide @ dump distarce that 1S small erougn to allcw the 
source tc Eke positiored many dump distances from the afex. 

Shaping the sand bottom into a ten degree wedge with 
a unifcrm and smootk interface proved to re atlceng and 
tedicus precess. ic facaditate Ehiscwenodeling, wooden 
supperts (two-by-fours) were mounted atony the length of 
froth sides cf the tark. These supports were elevated at cne 
end cf tke tank to achieve the required ten degree slope. A 
scraping device was ccnstructed with wooden supports anda 
metal scraping Elade that extended across the width cf the 
tank. This scraping device consisted of wooden sufforts 
along the top that reached across the tank and cculd ke 
Fulled aicng the elevated wocden supports on both side of 
the tank. This scrapirg device was pulled along the sufrorts 
repeatedly until a smooth slore of ten degrees was sculf- 
tured ircm the sand. 

Ecles had tc be drilled into the metal e¢craping 
tLlade Lecause when a solid Elade was used in the shallow 


Forticr cf the slope, water trap, ped behind the tlade was 
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Figure 3.2 Electronic Equipment Schematic. 


4e 


” 


forced remeath the Scraper gouging themwsaeer7 torre. ae 
Rottc i tcdeling techrigue waS Slow pecause, once a [ass was 
made cver the bettcn, the water became turctid and it was 
then it,cssible to see the LFotton. Hoen the bottom was rot 
visikle, it was imypcssible tc see where further suecthing 
waS necessary until the water settled several hours iater. 
Tae cit Jeu, as this smoothing process continued a silty 
residue Lecame separated frem the sand and settled cut on 
top cf the sand. Tkis residue would te easily resus;yended 
and eventually had tc fe removed using a water syphon. 


C. MEASTREMENT FROCELIORES 


Measuregent cf tke pressure field within the water sas 
done Ly licwering the receiver in depth at Specific ranyes of 
interest. The receiving hydrophcne was attacned to a rair of 
Ricrcweters at rigat angles to one another, that was in turn 
kolted tc a board which spanned the width of the tark. Once 
the Ecard was placed close to aranye of interest, che 
Micrcaeter was used to give fine adjustments in range anc 
the ctker in depth. 

The measurements were subject to both an accuracy anda 
Frecisicn error. On a given day, with the water level fixed 
and the cress-tank support set at a particular place in the 
tank, ait was possible to fosition the receiving hydrerkicne 
with an accuracy in depth and frange of plus or pinus 6.06 
centimeters {one turr of the micrometer). To prevent the 
sand ‘inland* of the apex frem drying out when not taking 
weasurements, enough water was added to the tank after each 
data run tckeer the sand ccrpletely submerged. The next 
time measurements were taken, water had to to ke reroved 
from the tank tc reestablish the beach. Because cr these 
small changes in the water level the hcrizontal fositicn cf 
the keack was subject tc a precison error estimated to be 


Withir plus O© Winws cue centinecare 
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Pre fital decisicn asscciated with the neasurererts 
Centeted ch whether tc use a triggered pulse cr a ccrtinucus 
wave (CR) Signal. With a triggered pulse it was possirle to 
distirguish the received signal from interference caused ry 
Moeeceicns Off the side of the tank. On the other harc, Ly 
using a triggered pulse there was a possibility that tne 
purse length was not long encugh for the acoustic energy 
meeocicteag with faths reflecting off the top and bettcm of 
Beem eater Column to cverlap the direct path from the scurce 
Mmemeeceiver. A CW signal would avoid octential pulse leaeth 
Beoplems, cut it would be impessible to Gistin,uish Letween 
tae actual signal ard interference. The use of a source with 
@ narrew hcerizontal ream reduced the erfects of reflecticns 
irom the side walls, kEut there was still the possibility o£ 
interference from side lobes reflected from the sides. 

It wéS Recessary to determine the best means tc take 
measurements. This sas done ky taking measurements at the 
same locaticn with different pulse lengths to determine the 
Ieguired pulse jlength to give consistent results and then 
Semparing these results to those obtained with a Ch saicnal. 
Toe tkird dump distarce (14.4 cm from the arfex) and jst 
fast the tenth dump Gistance (50.0 cm from the apex) were 
chosen, and measurements were taken for pulsed signals ci o4 
and z£6 cycles, andacC8 signal. These measuremerts can be 
Wennerigite 3.3 through Figure 3.5. In the figtrre 
Gepictirg resuits at the third dump distance, the Cu 
measurements are shown as the solid curve, the 64 txriygcer 
cycle results are disflayed as circular points, and the 256 
trigséer cycle results are depicted as triangular pecints. For 


clarity, the measurements taken just past the tenth dump 


Gistarce are shown ir two figures. Pagure 24.4 ccotpares Ch 
results (sclid curve) With the 64 triyger cycle sigral 
femeeular folnts). Figure 3.5 compares Ch ameasurenents 


(solic curve) with the 256 cycle results (circular [oints). 
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The pressure amplitude values were normalized ty dividing 
each [~resstre value ty the maximum pressure vaiue ir the 
itcmdwer die —Gepth values Were ‘also normalize@” bY GiVicing 
each certh Ly the maximum depth. 

The results in Figure 3.3 show gocd agreement Letweel. 
PieetbrEee curves. In particular, the Geseultsm@ecr ~the Ch 
Signal ard the 256 cycle signal are very sitilar. Ike 64 
cycle signal shows tke same behavior as the Ci sigral, fut 
GiIf£eEexrs 10 Lagnitude kelow rid-defth. 

The results at £C€.0 cm fircn the apex appear much nore 
Pomelacated than at the third dump distance. Figure 3.4 
compares results ackeived with the CW Sigcai anc the 34 
cycle sicnal. From the ficure it can be seen that the iwo 
sets cf results shcw good agreement. The measurerents 
depicted in Figure 3.£, which ccmpares the 256 cycie sigral 
and the Ch signal show foorer ayreenent, althcugh the 
general shapes cf both curves remain sizilar. 

Tkese results ixrdicate that interference from tne Side 
Peete SUtriclently smali sce that it 1S” possirie tc make 
measurezents with a Ck Signal to at least ten Gump distances 
from the apex. Long pulses could also te used, butetie 
difficulty cf making voltage reasurements with an cscillo- 
eeore CccCipared to™ reading a voltmeter, dictates that 


measurements should ke made with a CW signal. 
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IV. FODEL RESULT CCEEARISONS WITH LABORATORY Mraouktk 


A. ISTECLUCTION 


IfD model predictions and laboratory measuremerts were 
cbtaired fcr compariscn in EFoth a generai and detailed anai- 
ysis. The general analysis ccmpared results every five 
centineters from the beach to a range or 50.0 ca firci the 
keach. The measurements were spaced to give results at 
appreximately each cf the first ten dump distances. The 
detailed analysis cciryared results every centimeter freo 3.0 
to 17.0 cm from the ajyex. These measurements were taken to 
cbserve the sensitivity of the pressure field tc small 


changes ir range. 


Ee. TEE GCENERAL ANALYSIS 


Tke ceneral anlysis comrared IFD results asd lakcratcry 
measurements at apprcximately eacn of the first ten dump 
distarces (every five centimeters from the beach). tThese 
measurements were oktained by fixing the receiving hydro- 
Fhone sith respect tc the cress board and then moving the 
Foard out in increments of five centimeters to measure the 
desired field asa functicn cf depth at each range. the 
received CW Signal was read on the voltmeter. 

Tke Ccnfpariscns are shown in Figure 4.1 through figure 
4.10. In the Figures, the IFD predictions are disflayed asa 
solid line while the experimental measurements are Sacwn as 
circles ccnnected by a dashed line. Each pressure valué was 
nermalized to unity [Ly dividing it by the maximun [resstre 
value fcr the respective curve. The depth values were aiso 
normalized fy dividing each depth by the maximum depth in 


the water ccluon. 
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Compariscn cf£ Results at 1.0X. 
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Figure 4.6 Comparison of Results at 6.2X. 
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Figure 4.7 Ccmparison of Results at 7.3X. 
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Comparison of Results at 9.4X. 
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in general, the ;ressure patterns predicted ry the 1FD 
and tkcse measured, both become more compiicated as the 
momge Lrcmh the kteack increases. ea iensaal G os. there 1s 
Pimettative agreement in the Scale of the predicted features 
and tke scale of the measured features. Quantitative agzree- 
ment is lacking. This agreement in scale but not detail, 
Suggests that the phases amcng the normal modes predicted by 
the omacdel do not accurately  Leilect the, experimental 
erewaticr. 


Saab r EPETAILED ANAIYSIS 


Ike detailed analysis ccmpared IFD values wita ¢xcreri- 
mental measurements every centimeter from 3.9 cm to 11.0 ca 
from the Leach (0.7X tc 2.3X). fGhe laboratory neasuresents 
were taken by fixing the Foard at one Ilecsation azd tnen 
using the ricrometer to adjust the receiver to the desired 
range. Tke received Ch Signal was read from the voltmeter. 

ITbese results are depicted in Figures 4.11 tnrough 4.19. 
[iietne figures, the IFD predictions are shcwn as a so9lid 
line and tke experinental measurements are displayed as 
Circles. Each pressure value was normalized to unity by 
dividing it by the maximum fressure for the resrective 
curve. The depth values were normalized by dividing eaca 
depth ky the maximur cepth at the particular range. 

Tkere is gualitative agreement in the scale of the kLasic 
features fcr all rances, but yuantitative agreement is ret 


cbserved. The IFD patterns change very little thrceughcut the 


analysis, while the measured values chanyée more traricly 
(especially past 2.0%). Woe wesults “=atez- 1X and 253% 
(Figure 4.1&€ and Figure 4.19) Show that the pressure field 


Samco hance idlrly Significantly over a range as shott as cine 
centigzeter. At those ranyes where there iS poor agreement 


retween results, there is an Pidrecatm1On Or Fhase 
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Figure 4.11 Comparison of Results at 0.7X. 
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Figure 4.12 Comparison of Results at 0. 8X. 
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Comparison of Results at 1. 0X. 
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Figure 4.14 Compariscn of Results at 1. 3X. 
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Figure 4.15 Comparison of Results at 1. 5X. 
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Figure 4.17 Compariscn of Results at 1. 9X. 
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miter rcrence (Figure 4.18). FPOResd@scamecs 1655 tHaney sx, 
mmeary ELeCdicES that for a scurce at mid-aepth, the icwest 
Erorpagating mode shculd experience interference frem oniy 
the €vanescent tails of higher modes. Therefore, ir the 
megrcr Ci the detailed analysis, ore expects <« 
interference erfects. The experimental resuits in this 
regicn hecwever, shcw rather significant interference. fhe 
interference suggests unsuspected propagating modes are 
EEescrt, cc that tke evanescent tails are larger ttan 
empecteéd. from the trend in the curves it aprears as if the 
phase interference is not a tactor from 1.5X inward tcward 
the keack (Figures 4.11 through 4.15). 
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A. CCNCLUSZONS 
1. £eriormance cf the Jit 


the ITFD acoustic moiel did fot 


rh 


his andalys2s.0 
uncover any major failures cf model periormance 12 4 soe 
fied shallcw water environment. But this is not to say taat 
the mcdel consistently and accurately performs in such an 
cceanhn scenario. Althcuga comparison o£ IFD predicticns wita 
two ctker ocdels and with simple phySicat reasoning did pot 
uncover any inccnsisteéencies in performance, the agreement 
ketween IFC values and labcratory meaSurements is insutfi- 
client tc give complete confidence in tune rperforgrance cf 
€lther tre wodel or tke experizent. 

There 1S reasonable agreement between the scale of 
the features precicted py the IFD model and the experizert. 
Eoth the model and the measurements show increasing 
complexity as range is increased from the beach. But despite 
the sitilarities, tkrere 1S ECOr guaatitative agreemert in 
results. Cre possikle cause for the differences may re that 
the phases of the ncrmal mcdes predicted by the emcdel are 
extremely sensitive tc minor irregularities in the share of 
the interface and tke acoustic froperties of the Ecttcn. 
Frelininary work (LeSesne, 1984), suggests that the fhase 
relaticnshirs beéetweer modes is strongly dependent cn the 
distarce of the source from the apex even at great ranges. 
Consecuertly, it apjyears that the collective influence cf 
the ncernoal rodes 1s dependent upon careful geometric ccntrol 
cr the experiment. 

The detailed analysis arcund 2X, reveals that large 


phase interferences cccur where only one propagating mcde is 


expected. fas iItdicates Wiiemexicrenceomeof Gr extendsi 
SemcsCent tail withirt cutofivwes the Migher modes. Kavéesera 
ieee canneu, (1978), noted that the apparent phases c 
evanescent tails are extrenely sensitive to the 
MecmcrVitcRpnent, The resuits indicate that the ta 
decay guickiy, and its influence is pervasive 


and 4.3). From the trend in the curves, i 


ct 
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MieesemtntchiIcrence 12 hot a tactor any closer ¢ 
mmatienie > {figures 4.11 through 4.15). 


z- K&edelingsMeasurement Eroceiures 


Se SS ce — —_ = 


The iLabcratcry measurements represent an initial 
attermrt at modeling ar idealized shallow water 2nvircrment 
The €xperimental teckniy,jues are not without froblenms. All 
sets cfr measurements were repeated and snowed gcod a 
ment. The reproducibility cf the measurements suggest tha 
Tandcr é€rrers have [een minimized. However, systenratic 
Smears  Leliaty which also contribute to the dis 2 
cetween predicted and measured values. Some of these 
derive {£rcm ecuipment, such as the fact that the size cf the 
hydrcpfhone is of the same order of magnitude as the scale ce 
[eet iLicrouitimnne PLessure £Lields In addition, the envircn- 
ment in the tank may net be sufficiently close to the ideal 


envircnmgent assumed Ly the rodel. 


Fe. RECCHEBENDATIONS 
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is reccmmended. A cCcm@parison cf IFD predictions with 
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e 
models nct analyzed in this study, May uncover the cause 0 


the inccnsistency wken compared to measured values. A 
demaiied "study Or the phase interaction of the noraal ncd<cs, 
althcugh extremely ccaplex, say also offer insight intc the 


IFD perficriance. 
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It appears that the medeliny OF an idealizei Senator 
water Envirtcnment in the lakoratcry 15 O02 Varte 20 aeou 
modei verificaticn. L[espite the experimen rad2 tepewaao 
there waS gualitative agreement in tae basic features 
ketween aucdéel predicticns anad Jlakoratory weasureaients. Eut 


the experiment represents cnly an initial fro 


Intc He 


a) 


5 
rodelingymeasuremtent technigues of a sShailow water exvirc 
nent. Further refinerent of these tecnnigies and the =xieri- 
ental eéyguipment fay result 7G bettier my ene oe 
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FONNING THE CCNTOUR EFICT ON THE NPS COMPUTER 


Ae INTECLUCTION 


This Agpencix descriltes a fprocedure for running the MTL 
contcur flict on the NES computer. Detailed instructions for 


crunnirg the IFD ;rogram can be found in Jaéeger(1983). 


EF. CCEYING THE FILE FO USE 


Cnce tke TL values have FEFeen generated by the IFD 
prograr, aii that is needed to froduce a TII1 contour flict is 
the file FICIS FCRETRAN. This file is shcewn in Appendix E and 
can ke ccpied froaz a computer account maintained ry the 
Undersater Acoustics Curriculup. To link with this account 
and cktain a ccfy of the file, the user shculd freceed as 
fcllcws: 


(1) Meg7en tennineie 

(z) Enter: CP LINK) 0100e SIC tics. kie, 

(3) Wken prompted for the read passwcerd enter: UX . 

(4) Eften>s BEG igaee = 

(5) Enters COFY PiCiS FOatRANT CG —)—-aA = 
At this feint the PICIS FORTRANS file shouwic eres. cdomeecn tie 
user's A disk. 
C. RUNNING THE FROGEAM 


Eefcre runnirg tké€ program the user must obtain the TL 


values fircm the IFD frogram. These TL values must Le sert to 


pmecmauGaaekeiiecracr £6r the ELOTS program to be akle te tse 
Pelee lids Caneete Gdc¢re by  ~-lecing a WRITE Statement in the 
IFD p~regram that sends the values to a data disk. The sgodi- 
fied JFD déficted in Appendix A uses this technigue anc can 
ke used as a guide. 

Croce tke data disk has teen created, the user =uotrst 
assicn téempcrary disk space (TIDISK) to give the precgran 
sufficiert room to generate the IL contcurs. For mere infcr- 
Maticn cnr how to assign temporary disk space, s€¢e IAPS 
Technical Ncte TN-VM-C1 which is availible in the ccomruter 
eomsultart's office. 

With these initial steps ccmpieted, all that remains is 
to cetrile and run tke program. The program can be ccnrziled 


ey: 
EMters Fr CHI Ci VFLOTS . 


The Eregraz must be run at the TEK618 graphics terminal 


Updenr  DaSst lA. §inis Can be dcne Ly; 
EntCr ee loorLA < 


The uséxr will then Le precmpted for the compiled fFertran 
Peegreterame sdnd the rile definitions for the data cisk, 


kefore the program wiji run. 


SOUBKBCE LEPTIH SENSITIVITY ANALYSIS 


During the laboratcry experiment, two sueniemeniee sets 
cf mé€asurements were taken to obtain an indicaticr cf the 
Sensitivity of the pressure auplitude to changes in the 
source depth. The first set of measurements was okEtained by 
varying the source depth swith the receiver fixed at 1.0X 
(4.8 cm fren the beach), and lowered in depth to the Ecttcn. 
For tke seccnd set of meaSurerents, The recéiver was fixed 
at the third dump distance (14.4 cm from the beach), and 
measurements were taken with the source fixed at 5, 7, 9, 
11, 12, 15, 17, 19, 2, 23, Ze; 21,0 277 end and Sacer mee 
the surface. 

The results of the first analysis are shown in Figtre 
De1. In the figure, fressure was normalized by dividing by 
the nwaxigtum pressure and depth was normalized by dividirg by 
the depth cf the water column. Although theory predicts only 
che proyagating mode at this distance, the results show 
modal phase interferelce. 

The results cf tke seccnd analysis are shown in Figtre 
D.2 tkrcough Fagure D.€. In the f£qures, “cach Curve 7e7ce 
sents a set of measurements taken with the source fixed at a 
specific defth. For convenience, more than one curve is 
shown cn a sgiven figtre. All fressure values were norgalized 
Ey dividing by the tke maximur fressure in the field. the 
depth values were necrmalized Fy dividing by the maxiuzun 
water derth. 

lregeneral, at =:.0X there appears to be modal fhase 
interference which is influenced Ly source derfth. Giver the 
rougkness of this analysis, the details of this interference 


are okscure. However, the modal interference does net afpeéar 
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to Le incorsistent sith either the placement of the scurce 


cr with tne previous exfperinaental measurements. 


ile 
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Figure wD. 1 source Sensitivity Analysis at 1.0X. 
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Figure D.4 Measurezents with Source Depth of 17, 19, and 21 Cn. 
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